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Abstract 

We study conditions under which a partial differential operator of arbitrary order 
n in two variables or ordinary linear differential operator admits a factorization with a 
first-order factor on the left. The factorization process consists of solving, recursively, 
systems of linear equations, subject to certain differential compatibility conditions. In 
the generic case of partial differential operators one does not have to solve a differen- 
tial equation. In special degenerate cases, such as ordinary differential, the problem is 
finally reduced to the solution of some Riccati equation(s). The conditions of factor- 
ization are given explicitly for second- and, and an outline is given for the higher-order 
case. 

1 Introduction 

Factorization of linear differential operators (LDO) is a very well-studied problem and 
a lot of pure existence theorems has been proved (e.g. pQ, |Hj, [HI)- The case of lin- 
ear ordinary differential operators (LODO) has been studied in more detail; various 
algorithms are known for factoring LODOs over different differential fields, e.g. [3], 
jllj and others. Probably the first algorithm for the case of the simplest possible 
differential field (rational functions) is described in [2]. In this case factorization into 
irreducible factors is known to be unique [Jj and in general, the problem finally 
reduces to the solution of a Riccati equation. 

Much less is known about the factorization of linear partial differential operators 
(LPDOs) which would appear to be much more complicated than LODO. In several 
variables the naive definition of factorization as representing a given n-th order operator 
as a composition of lower order operators does not have a nice uniqueness property; 
this is illustrated by an example of E. Landau (formulated for left composition) 

(d x + i)(d x + i)(d x + xd y ) = [d 2 x + xd x d y + d x + {2 + x)d y ](d x + 1). 

Nevertheless, some interesting results have been obtained also in this direction; re- 
cently Grigoriev and Schwarz 4 have given an algorithm for factoring an LPDO with 
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separable symbol. Some generalized definitions of factorization gave rise to another fac- 
torization algorithms for factoring LPDO and for factoring of systems of LPDEs 
with rational function coefficients ^3], |14) . 

It seems to be a general opinion that the naive factorization - even for second- 
and third order differential operators in two variables - requires "the solution of a 
partial Riccati equation, which in turn requires to solve a general first-order ODE 
and possibly an ordinary Riccati equation. The bottleneck for designing a factoriza- 
tion algorithm for LPDO is the general first-order ODE which makes the full problem 
intractable at present because in general there are no solution algorithms available" £Q . 

We consider it possible to construct an explicit algorithm which could be used for 
absolute factorization of LPDO of arbitrary order n in two variables into linear factors. 
The word " absolute" means that we do not fix the coefficient field from the very be- 
ginning and that our only demand on the coefficients is that they be smooth, i.e. they 
belong to an appropriate differential field. The procedure proposed here is to find a 
first order left factor (when possible) in contrast to the use of right factorization, which 
is common in the papers of last few decades. Of course the existence of a certain right 
factor of a LPDO is equivalent to the existence of a corresponding left factor of the 
transpose of that operator, so in principle nothing is lost by considering left factoriza- 
tion. Moreover taking transposes is trivial algebraically, so there is also nothing lost 
from the point of view of algorithmic computation. 

The use of left factorization makes all the algebraic calculations needed much eas- 
ier. It also yields explicit necessary and sufficient conditions on the coefficients of the 
operator for the existence of a factor. The simple form of these conditions allows us 
one to see the principal difference between factorization of LPDOs and LODOs - in 
the first case the factorization problem can be solved by pure algebraic methods in the 
"generic" case; the LODO is one example of the various degenerate cases in which it 
is necessary to solve a Riccati equation in order to factor in two variables. 

In Sec. 2 the general idea of factorization (in the generic case) by purely algebraic 
methods is presented and explicit conditions of factorization for an LPDO of order 2 
are written out. It is also shown how the possibility of pure algebraic factorization for 
LPDO does not contradict the necessity to solve a Riccati equation in order to factor 
an LOPO. Similar results for LPDO and LODO of order three are demonstrated in 
Sec. 3, while in Sec. 4 the general procedure for factorization for LPDO of arbitrary 
order n is presented. Some interesting examples are constructed in Sec. 5 and Sec. 6 
contains a brief concluding discussion. 



2 LPDO of order 2 in two variables 



We consider an operator 





with smooth coefficients and look for a factorization 



M = {pid x + p 2 d y + pz){pid x + p 5 d y + p e ). 
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Let us write down the equations on pi explicitly, keeping in mind the rule of left 
composition, i.e. that d x (ad y ) = d x (a)d y + ad xy . 

Then in all cases 

O20 = PlPA 

Oil = P2PA + PlPb 

O 2 = V2Vo 

aio = £{p±) +P3P4 +P1P6 

OOl = £(Pb) +P3P5 +P2P6 

^a 00 = C(p 6 ) +P3Pe 
where we use the notation C = pid x + P2d y . 

We begin with some preliminary remarks. To some extent the factorization problem 
is to be treated as a local problem; for example if some coefficient does not vanish 
identically, we may wish to restrict to a region where it has no zeros. Of course 
insofar as the formulas obtained are analytic functions of the data, if the data are also 
analytic then the formulas may continue globally. With this in mind, we note that 
if the operator is not (globally) of order one, then we may, after a linear change of 
variables if necessary, assume that 

020 ¥= 0. 

Then necessarily pi ^ 0, and we may assume without loss of generality that 

Pi = 1- 

Then the first three equations of 2SysP, describing the highest order terms are equa- 
tions in the variables P2,Pa,P5 and to find them we have to find roots of a quadratic 
polynomial. Having computed P2,Pa,P5 one can plug them into two next equations of 
2SysP and get a linear system of equations in two variables P3,P6 which can easily be 
solved if a; is a simple root. Then the last equation of 2SysP will give us the condition 
of factorization. 

In fact, at the first step from 

O20 = Pi 

< an = P2P4 + P5 
, a 2 = P2P5 

it follows that 

V 2 (-P2) ■= a 2 o{-p2) 2 + au(-p 2 ) + a 2 = 

The choice of different roots of the characteristic polynomial V2 gives us different 
(potential) factorizations of the initial operator. (As we shall see, there may be no 
such factorization, or there may be one with one choice of root and not with another 
choice.) Let uj be a root of V2, and take 

p 2 = -uj, C = d x - tody. 

This leads to a linear system for P4, p^ with uj as parameter: 

Pa = a 2 o, 
-ujpi +p 5 = a u 



{2SysP) 
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or, in matrix form, 



" 1 0" 




Pi 




O20 




Pi 




1 0" 




O20 


-L0 1 




P5. 






1 


P5. 




u 1 




an 



Thus 



pi = 1 

P2 = ~U 

Pa = a 2 o 

p 5 = a 20 u; + an 



(2Pol) 



At the second step, from 



\a w = £(p 4 ) +J> 3 P4 +P1P6 
\a 01 = £(p 5 ) +P3P 5 +P2P& 

and from (2 Pot) we get 

jaio = £a 2 +P3«20 + P% /or - \ 

< (2Lin) 
[a i = £(an + a 2 o^) +P3(an + O20^) - ^P6- 

If the root u) is simple, i.e. V' 2 (oj) = 2a 2 $u) + a n / 0, then these equations have the 
unique solution 

Loaio + aoi - w£a2o — £(a 2 ow + ai1 ) 
P3 = 



P6 



2a2o^ + an 

(a 2 ouJ + an)(aio - £)a 2 o - 020(^01 - £(a20^ + «ii) 



2a 2 o^ + an 

At this point all coefficients pi,p 2 , ■■■,Pe have been computed. 
At the third step from 

a 00 = £{pe) +P3P6 
the corresponding condition of factorization can be written out explicitly: 



a _ £ [ uai ® + a 0! ~ ^( 2q 20^ + on) ) | ^aio + Qqi ~ £(2a 2 o^ + an) ^ 
\ 2a 2 o^ + an J 2a 2 o<^ + an 



ZU^ I "11 J ^i*2U w I "11 

a2o(Qoi ~ £)(«20^ + Qll) + (fl20^ + Qii)(aio - £a2o) 
2a 2 o^ + «ii 



and we have just to check it. Thus factorization in two variables, when possible, can 
(generally) be accomplished by purely algebraic means. 

Remark. In the preceding derivation we assumed that to is a simple root so that 
(an + 2a2o^) / 0. Suppose now that the condition fails identically in some region. 
Then by eliminating p$ from (2Lin) we obtain the necessary condition 

aiokJ + aoi — C{a 2 QUj + an) — u}£a 2 o = 0. (2Degen) 
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If this condition is satisfied, then we may use the first equation of (2 Lin) to determine 
Pq from p3. Then the last equation of (2SysP) gives a Riccati equation for p%. 



Let us look explicitly at the case of an ordinary differential operator of second order 
as a particular case of our initial LPDO. Here to is a double root of the characteristic 
equation and the necessary condtion (2Degen) is necessarily satisfied. In fact 2SysP 
takes the form 

'020 = Pa 

010 = £(P4) +P3P4 + P6 (2SysO) 
a o = C(p e ) +P3P6 

i.e. we have only three non-trivial equations. The first equation shows again that P4 = 
020 ! £-Pa = fl 20 is known, and p§ can written out explicitly as p% = (aio — 020 — 02o)P3- 
Substitution of this expression for p^ into the equation for aoo gives us immediately 

OOO = OlO - O20 - O20P3 - O20P3 + 0l0 ~ ^20 ~ 020^3)^3 (2) 

which is the Riccati equation 

; , / 2 , 2ci2o — oio , . °oo + «20 — «io n r , 
ip + ip H V "I = U, tor P3 = ip. 

«20 ^20 

In order to be able to find the general solution explicitly one has to know at least one 
of its solutions. Of course the solution will depend on the explicit form of = aij{x). 



3 LPDO of order 3 

Now we consider an operator 

A 3 = ^ a jkd J x dy = a ZQ dl + a 2X d 2 x d y + a 12 d x d 2 y + a 03 dy 3 

j+k<3 

+a 2 odl + aud x d y + a 02 dy + a w d x + a 01 d y + a 00 . 
with smooth coefficients and look for a factorization 

A3 = (pid x + p 2 d y + pz){pAdl + p^d x d y + p 6 dy + p 7 d x + p 8 d y + p 9 ). 
The conditions of factorization are described by the following system: 



O30 


= PlP4 


021 


= P2P4 + PlP5 


Ol2 


= P2P5 + PlP6 


O03 


= P2P6 


O20 


= £(p 4 ) +P3P4 +P1P7 


Oil 


= £(p 5 ) + J> 3 £>5 + P2P7 + PlPS 


O02 


= £{P6) +P3P6 +P2P8 


OlO 


= £{P7) +P3P7 +P1P9 


aoi 


= C(p 8 ) +P3P8 +P2P9 


000 


= jC(p 9 ) +P3P9 



(3) 



(55ysP) 
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with C = pid x + P2d y . 



Once again we may assume without loss of generality that the coefficient of the 
term of highest order in d x does not vanish, and that the linear factor is normalized: 

030 + 0, pi = 1. 

The first four equations of SSys Var describing the highest order terms are equations in 
the variables P2,Pa,P5,P6- Solving these equations requires the choice of a root — p2 of 
a certain polynomial of third degree. Once this choice has been made, the remaining 
top order coefficients Pa,P5,Pg are easily found. The top order coefficients can now 
be plugged into the next four equations of SSysP. The first three of these four equa- 
tions will now be a linear system of equations in the variables P3,P7,ps which is easily 
solved. The next equation is now a linear equation on variable pg which means that 
all variables pi,i = 1, 9 have been found. The last two equations of SSysP will give 
us then the conditions of factorization. 

Namely, at the first step from 



O30 
021 

ai2 
ao3 



PA 

P2PA + P5 
P2P5 + P6 
P2P6 



it follows that 



Vz{-P2) ■= a 3 o(-p 2 ) 3 + a 2 i(-p 2 ) 2 + a 12 (-p2) + a 3 = 0. 

As for the case of second order, taking p2 = -co, where a; is a root of the characteristic 
polynomial we get a linear system in P4,P5,Pe with to as parameter. Then again 



P2 = -oj, C = d x - Lvdy, 



which leads to 



O30 
1,012 



PA 



-Ujp 4 + p 5 

-up 5 + p e 



i.e. 



1 

-oj 1 
-oj 





P4 




030 




pa 




' 1 





0" 




030 




P5 




o 2 i 




P5 




OJ 


1 







021 




Pa. 




Oi2_ 




P6. 






OJ 


1_ 




CL\2_ 



Thus 



Pi 

P2 

Pa 
Pb 

{P6 



1 

—oj 

O30 

a 30 oj + a 2 \ 

a^ooj 2 + 021^ + a\2- 



At the second step, from 



(3PoZ) 
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'a 20 = £(p 4 ) 

< an = £(p 5 ) + P1P& 

^02 = £ (p 6 ) + P3PQ + P2P8 

and (2 Pot) we get 



a 2 o - £a30 = P3O30+P7 

an - C(a 30 uj + a 2 i) = p 3 (a 30 uj + 021) - up 7 + p 8 (3Lin*) 

a Q2 - C(a 30 uj 2 + a 2 iw + ai 2 ) = p 3 (a 30 u 2 + a 2i a; + an) - ojp 8 . 

As a linear system for p 3 , P7, ps this has determinant 

3a 30 w 2 + 2a 2 i^ + ai 2 = V(u>), 

so if uj is a simple root the system has unique solution 

w 2 (a 20 - £a 30 ) + u(au - C(a 30 uj + a 2 i + a 02 )) - £(a 30 w 2 + a 2 \uJ + a i2 ) 



£»3 



3a 30 uj 2 + 2a 2 iu; + a i2 

_ Q20 — £^30 ^30 

3a 30 u; 2 + 2a2iw + an 3a 3 ow 2 + 2a2iw + an 



uj(a 2 o - £a 3 p) + an - £(a 3 pu; + a 2 i) a 3 puj + agi 

3a 30 w 2 + 2a 2 iw + a i2 3a 30 ^ 2 + 2a 2 i^ + a i2 



In order to find the last coefficient pg we use the next equation of (3SysVar), namely: 

aio = C(p 7 ) + p 3 P7 + P1P9, (3Lin * *). 

At this point all coefficients pi,i = 1, ...9 have been computed, under the assump- 
tion that uj is a simple root. 

At the third step from 



a 01 = C(p 8 ) +P3PS+P2P9 
a o = C(pq) +p 3 p 9 

all the necessary conditions for factorization can be written out. We do not do so 
here because the formulas are tedious and do not add anything to understanding the 
main idea. If the conditions are satisfied, the explicit factorization formulae could be 
written out as for the second-order operator. The difference is that in this case the 
polynomial defined by the highest order terms is of degree 3 and we have not one but 
two conditions of factorization. 

Regarding now LODO as a particular case of LPDO let rewrite SSysP with p\ = \ 

as 

030 = Pa 

O20 = P3P4+P1P7+P4 
aiO = P3P7 + P9 + P7 
^00 = P3P9+P9- 



(SSysO) 



7 



As in the second order case oj = is a multiple root of Vz . After eliminating P4 = 030 
and solving the second and third equations for pi and pg, respectively, in terms of P3, 
we see finally that the factorization problem for a third-order LODO is equivalent to 
the following system of first order nonlinear ODEs for the single unknown function ^3 : 

aio = P7 (a-20 ~ a-30 ~ Pi) 1 030 + P% + Pr 

a-00 = P9 (O20 - ^30 - Pi) I O30 + P9- 

or, equivalently, a certain second-order nonlinear ODE. 



4 LPDO of arbitrary order 

Now consider the factorization problem for the LPDO of order n in the general form 
A n = Y, a jt ^ = (pA+M+P3)( I'j^i) ( 4 ) 

j+k<n j+k<n 

and use the same assumptions and notations as above: 

• some ajk with j + k = n is not = 0, so 

• after a linear change of variables one can assume that a„o 7^ 0, 

• set p\ = 1 and pjk = if j < or k < 0, 

• denote P2 as —to and C = p±d x + P2d y = d x — ojd y . 

Then the equations describing terms of highest order take the form p±Pj-i ) k+P2Pj,k-i = 
ajk, i-e. 

Pj-i,k ~ upj,k-l = ajk, j + k = n, (5) 
while the equations of lower order take the form 

£pjk + PzPjk + Pj-i,k ~ MPj,k-i = a jk, j + k <n (6) 
Keeping in mind that j = n — k for the first n equations (jSJ let us rewrite them as 

Pn-l-k,k — <^Pn-k,k-l = a n-k,k, fe = 0, 1, n 
or in matrix form as 



1 








... 0" 




Pn-1,0 




a n ,o 


—to 


1 





... 




Pn-2,1 









—uo 


1 


... 




Pn-3,2 




a n -2,2 










-to 1 




jP0,n-l_ 




_ a l,n-l_ 



If we knew u, the unique solution could be written out as 



Pn-1,0 




1 








... 0" 




fln,0 


Pn-2,1 




U! 


1 





... 






Pn-3,2 




cu 2 


a; 


1 


... 




a n -2,2 


jP0,n-l_ 




0J n ~ l 


oo n ~ 2 




to 1 




_Gl,7l-l_ 
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so 



k k—1 

Pn-k,k = a n ou) + a n -i,iw H + a n -k±- 



(8) 



The equation with k = n, together with u;po,n-i = «o,n gives us finally 

V n (v) := a n ,ou n + an-i^w" -1 + ... + ai >n _iw + a , n = 0. 

Choosing a; to be any root of this polynomial, we get then a unique determination of 
Pj t k for j + k = n. 

Consider now the n equations of (JHJ) corresponding to the terms of order n — 1, i.e. 
the case j + k = n — 1 : 



Pn-2-k,k - Wp n -k-l,k-l = b n -l-k,k — PzPn-X-k,ki 

where the functions p n -\-k,k an d 

b n -l-k,k = a n-l-k,k — £Pn-l-k,k 

are known already. In matrix form these equations are 



0) 



Pn-1,0 


1 





... 


" 




P3 




frn-1,0 


Pn-2,1 


— U> 


1 


... 







Pn-2,0 




bn-2,1 


Pn-3,2 





—UJ 


1 ... 







Pn-3,1 




bn-3,2 


_P0,n-l 








... 


—to 




_P0,n-2_ 




pQ, n -l 



Inductively, the determinant of a matrix of this form is 

(-I)"' 1 [pn-l^- 1 + + • • • + P0,n-l] 

Equations © imply that the term in brackets is 



ndnfiU 



n-l 



+ (n - l)a n _ 2 ,iw n - 2 + • • • + ai,„_i = ^(w). 



Therefore if a; is a simple root, the equations © for ^3 and the p n _2-k,k have a 
unique solution. Multiplying the j-th equation of the system by ui n ~ J and adding, we 
determine 



V' n {oj)p ?J = bn^flOJ 11 ' 1 + b n -2,lOJ n ~ 2 + ■■■ + b , n - 



(10) 



If a; is a simple root, then this equation determines ps, and we may determine the 
Pni2-k,k by rewriting the first n — 1 of these equations with the p% terms on the right 
hand side, so that the solution for the p n -2-k,k m terms of p% and other known functions 
is 



Pn-2,0 




1 





Pn-3,1 




UJ 


1 


Pn-4,2 




OJ 2 


u 


J>0,n-2_ 






u n-3 



... 

... 

... 

uj 1 



bn-1,0 ~ P3Pn-l,0 
bn-2,1 ~ P3Pn-2,l 
b n -3,2 ~ P3Pn-3,2 

h,n-2 ~ P3Pl,n-2. 



(11) 



as before. 



Continuing to assume that uj is a simple root, let us follow this process one more 
step. The first n — 2 of the n — 1 equations for the coefficients p n -3-k,k can be written 
in a matrix form like (|7j). with known functions on the right hand side, so the solution 
can be written down as above. Thus at this point we have determined the coefficients 
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Pj of the first factor and all coefficients pjk, j + k < n — 3 of the second factor. However 
there is one more equation at this step: 

-Wp ,n-3 = °l,n-3 ~ (£ + Pz)Pl,n-3- 

Since both sides have already been determined, this is a necessary condition for the 
existence of a factorization with this choice of root u. Continuing, the same situation 
occurs at each step. Thus if the condition is satisfied at one step we may proceed to the 
next step, and either the process fails or there are exactly n — 1 polynomial conditions 
that are satisfied by the coefficients, the (simple) root u, and their derivatives. 

Suppose we drop the assumption that u> is a simple root and assume that V n {oj) = 
in some region. Then the equation with zero right hand side, is the necessary 

and sufficient condition for solvability of the full system of equations for p^ and the 
p n -2-k.k- If this condition is satisfied, then we may make any choice of p% and de- 
termine the p n _2-k,k from 1)11(1 . At subsequent stages we will encounter factorization 
conditions that are expressible as polynomial equations in p% and its derivatives (with 
coefficients that are polynomials in the ajk and their derivatives and in the root u and 
its derivatives); these can be thought of as generalized Riccati equations for ^3. 

In summary: starting with a simple root of the polynomial V n {io) in the generic 
case a n fl 7^ 0, the factorization process proceeds algebraically (in the coefficients cijk, 
the root u>, and their derivatives) as far as it can, as determined by n — 2 polynomial 
obstructions. Starting with a double root, then subject to one such polynomial con- 
straint, the remaining constraints are n — 2 generalized Riccati equations for p^. 

We call this eliminating procedure order-reduction for obvious reasons. In trying 
to factor an operator of order n as a product of a first order operator and an operator 
of order n — 1 , we look at the equations associated to the terms of given degree of the 
operator of order n, in descending order. The first set of equations gives us the top 
order terms of both factors (modulo choice of a root of the characteristic equation), 
the next set determines the next lower order coefficients, and so on. 

5 Some examples 

5.1 General second-order LPDO 

• Example 5.1.1. Hyperbolic operator: 

020 = 0, an = 1, 002 = 0, i.e. L 2 ■= d x d y + a w d x + a m d y + a o- 
In this case either of two different factorizations may be possible: 

- (5.1.1.1 ) if a o = d x a w + a w a i, then 

L 2 = (d x + a i) (d y + aio) 

oi- 



lO 



— (5.1.1.2 ) if aoo = d y ao\ + aioaoi, then 

L 2 = (d y + aio) (d x + 001) • 
Consider the case of such an operator known as a Poisson operator: 

L 2 = d x d y H — (ad x + (3d y ) + - — ^— ? 

where a,/?, 7 are constants. 

Due to (5.1.1.1) and (5.1.1.2) factorization is only possible if 7 = a(/3 — 1) or 

7 = /3(a — 1) and its explicit forms are, respectively: 

L 2 =(d x + ^-)(dy+ " 



x + V J V x + y 
or 



x + y J \ x + y 

We note that in the particular case 7 = 0, the Poisson operator reduces to the 
operator considered in Example 4 from [1] and the results there can be obtained 
immediately from the formulas above. 

• Example 5.1.2. Parabolic operator 

The operator 

A p = a 2 odl + a u d x dy + a 02 dy + a w d x + a 0l d y + a 00 

is parabolic only if a\ 1 — 4a2o«o2 = 0. This condition means that the characteristic 
polynomial 

V 2 := a 2 QUJ 2 + anuj + 002 

has a double root u = —an/a 2 Q and it follows from the discussion in Sec. 2 that 
factorization requires solution of a Riccati equation for ^3. 

• Example 5.1.3. Elliptic operator 

The operator 

A e = a 2Q d 2 x + a n d x d y + a 02 dy + a w d x + a id y + a 00 

(with real coefficients) is elliptic by definition if a\ 1 — 4a2oao2 < 0, so that the roots 
of the characteristic polynomial V 2 are complex. Therefore, the factorization is 
only possible over a field that contains C. 

5.2 Role of zero term; order of factors 

Consider the following second order operators in two variables: 

A a = d 2 x -d 2 y + xd y + yd x + \(y 2 - x 2 ) + a, (12) 
where a is a real constant. 
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• Example 5.2.1 

We look for a factorization in the form 

A-a = (pid x + p 2 d y + pz){pid x + p 5 d y + pq). 

Necessarily p\Pi = 1, and with no essential loss of generality we take p\ = p^ = 1. 
Then it is not difficult to show that the factorization can only have one of two 
forms: 

A a = L+L~ = (d x + d y + b)(d x -d y + c) (13) 

or 

A a =L~Lt = (d X -dy + b)(d X +dy + c) (14) 

It follows from the results above, and can be checked directly as an exercise, that 
the following is true: 

— The operator 

A a = d 2 x - dy + yd x + xd y + \{y 2 - x 2 ) + a, 

does not have a factorization of either form Q13JI or (|14|) unless a = dbl. 

— The operator 

A x = d 2 x - d 2 + xd y + yd x + \{y 2 -x 2 ) + l 

has a factorization of the form ()13|) but has no factorization of the form (I14|) . 

— The operator 

A_a = d 2 x - d 2 y + xdy + yd x + \{y 2 - x 2 ) - 1 
has a factorization of the form (|14|) but has no factorization of the form (|13[1 . 

In fact 

A x =[d x + d y + \{y- x)] [d x -d y + \{y + x)] ; (15) 

and 

A-i = [d x -d y + \{y + x)] [d x + d y + \{y - x)] . 

This example shows that the existence of a factorization depends crucially on the 
term of order zero, and also that an operator may have a left factor of a certain 
form but no right factor of the same form, or conversely. 

As we have shown, in two variables the question of existence of a first order left 
factor of a given form can be settled, and the the factors calculated, in a system- 
atic way using procedures that are simple and (in general) purely algebraic. 

• Example 5.2.2 

The question of existence of a first order right factor can be reduced to the 
question for a left factor by taking the formal transpose. As an illustration, let 
us find a right factor for operator A\ in this way. In order to take advantage of 
what is already known, we make a coordinate change y — ► —y, which converts A± 
to 

Al =d 2 x -8 2 - yd x - xdy + \(y 2 - x 2 ) + 1. 
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The transpose is 

A\ = d 2 x - 8 2 + yd x + xd y + \(y 2 - x 2 ) + 1 = Ay = L+(L^). 

Therefore 

A = (LD*(L+)' = (-L^n-L+f. (16) 
Now according to lfT3|) . 

(-L+Y = - [d x + 8 y + \{y -x)] t = d x + d y - \{y - x). 

Undoing the coordinate transformation y — > —y converts this last operator to 

d x ~ d y + \{y + x). 

According to (|16|) . this is a right factor for Ay, which is consistent with Q15|). 

5.3 Classes of factorizable LPDOs 

Putting some restrictions on the Ojj as functions of x and y, one can describe all 
factorizable operators in a given class of functions. Let us exhibit explicit factorizations 
for one case of hyperbolic operators: 

A := d xx - d yy + a w d x + a idy + a 00 , (17) 

i.e. a 2 o = l,an = 0,a 2 = -1. 

In this case there are only two possible factorizations: 

• case 1: if 4aoo = 2(d x + d y )(a±o + aoi) + (af — Ooi) holds, then we have the 
following factorization: 

a /a i fl i fl W" a 0l\ /o a i a 10 + a 0l\ / 10 x 

A 2 =\d x + d y + J\d x -dy + I (18) 



• case 2: if 4aoo = 2{d x + d y ){a\Q — aoi) + (af — «oi) holds, then we have the 
following factorization: 

a (a , a i ai + fa a , a 10~ a 0l\ /ln x 
A 2 =l^ + ^ + II I. (19) 



Let us consider as a simple example the case when the lower order coefficients 
aio, aoi) «oo are linear functions of one variable x, i.e. 

oio = »io,i^ + »io,o, 
< aoi = &oi,i» + 6oi,0) 

Ooo = ^00,1^ + ^00,0- 

Then the condition of factorization (24) is 

46oo,ix+46 o,o = 2b 1 o tl + 2b i :1 +bl 0tl x 2 + 2b 1 o :1 b 1 o fi x + bj 0fi -bl 11 x 2 -2boi,ib olfi x-bl lfi 
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which leads to 



zero term: 



linear term: 



quadratic term: 



4&oo,o — 2&io,i + 2&oi,i + bf - 6qi n, 



4&oo,i — 26io,i&io,o — 26qi, lftoi.o^ 



= &? 0) i " &oi,i- 



The last equality gives immediately 610,1 = ±&oi,i- Let us consider, for instance, the 
case &io i = — boi i = b. Then 



4600.0 - ^10,0 ~ b oi,o 

26 0.1 = -&0io,o + bo^o), 



or in parametric form: 



601.0 = V%h sinhi 2 , 
ho,o = V%h coshi 2 , 
froo,o = \t\i 

600.1 = -\/2tit 3 (sinht2 + cosht 2 ), 

&10,1 = *3 
&01,1 = — *3 

with arbitrary real parameters t±, t 2 , t% G -R. Varying these parameters one can get coef- 
ficients &oi,0) ^io,0) ■•■) ^00,0 and, therefore, coefficients of the initial operator a±o, aoi, aoo, 
that belong to some prescribed field. For instance, with ti,t$ € Q and i 2 = 0, 



601.0 = 0, 
bio,o = V%h, 
boo,o = \t\, 

600.1 = — V^tlt3> 
&10,1 = *3 

&01,1 = — *3 



and we get 2-parameter factorization in Q{y/2): 
A 2 = \d x + d y + 



d x -d y + 



V2h 



The case &io,o = &oi,o can be treated in a way analogous to the previous one and case 2 
in a way analogous to case 1. Combining all the results one will give a description of all 
factorizable LPDO of the form (22) with corresponding restrictions on the coefficients. 
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6 Conclusions 



The elimination order-reduction procedure has been presented for explicit factoriza- 
tion of LPDO of arbitrary order. It was shown that in generic case, i.e. in case 
when characteristic polynomial has a distinct root, this procedure is wholly algebraic 
and no differential equation has to be solved. Otherwise the factorization problem is 
equivalent to solving some Riccati equation. It was also demonstrated that different 
types of LPDO have different factorization properties, namely, algebraic factorization 
of hyperbolic operators, when possible, could be done over 1Z, while factorization of 
elliptic operators is only possible over C. In case of parabolic operator no algebraic 
factorization is possible. 
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